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1. Introduction 

There is a renewal interest in two, in principle, different topics. One is the multi- 
symplectic formalism for classical field theories, introduced in the latest seventies by the 
Polish school led by Tulczyjew (see 0] for details), and independently, by P.L. Garcia et al 
H |l| , and H. Goldsmichdt and Sh. Sternberg @ (see §, H IDf 111' 111' U H' HI for 
recent approaches). The other topic is mechanical systems subjected to non-holonomic 
constraints. This last topic, being a classical subject in the literature, has recently be- 
come very prominent in the area of Geometric Mechanics. 

Our aim in this paper is to give a multisymplectic setting for first-order field theories 
subjected to non-holonomic constraints, that is, the lagrangian depends on the first 
derivatives of the fields, and, in addition, not all the values are allowable. We do not 
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know interesting examples in physics: for instance, div E = in electro-magnetism, 
or div E = p in clcctromagnetism interacting with charged matter (see [p7| ) are non- 
holonomic constraints, but they are dinamically invariant, so that nothing new can be 
added using our formalism. Continuum Mechanics seems to be the main and most 
interesting source of examples. 

In H^] the authors consider incompressibility constraint in fluids and elastic solids, 
but this type of constraints become holonomic after integration in the space of Cauchy 
data. A genuine type of constraints are exhibited by an elastic body rolling with constant 
velocity, while deforming (for instance, an automobile tire on pavement ]2§|]). Another 
important family of constraints are those appearing in the dynamics of media with mi- 
crostructure. Let us recall that a simple body is mathematically modelled as a three- 
dimensional manifold B which can be covered by just one chart; therefore, a configura- 
tion is just an embedding of B into M 3 , and a deformation is a change of configuration. 
A medium with microstructure is more involved: it is modelled as a principal bundle 
P — ► B which is embedded in a physical ambient principal bundle. The configurations 
and the deformations are now principal bundle isomorphisms Q . A multisymplectic set- 
ting for elastic media was developed in p2] | , and it is a natural continuation of the nice 
work in p6[ ] . Many of these continuum systems exhibit kinematic constraints given by 
the existence of a connection in the ambient principal bundle |3| ^) . 

In this paper, apart from presenting a multisymplectic setting for first-order field 
theories subjected to non-holonomic constraints, we also develop the corresponding infi- 
nite dimensional setting in the space of Cauchy data, proving that the constraint forces 
can be integrated over the Cauchy surface to get a constrained 2-form on the space of 
Cauchy data. So, the equations of motion are obtained modifying the usual 2-form with 
this constrained 2-form. 



2. Multisymplectic geometry 

In this section we will recall some fundamentals on multisymplectic manifolds which 
will be needed later. We refer to |J for more details. 

Definition 1. A closed fc-form cu on a manifold M is called multisymplectic if the 
mapping 

X e T X M ^i x tu e (A fe_1 .E% 

is injective for all x G M. In this case, the pair (M, ui) is called a multisymplectic 
manifold. 

As is well-known, cotangent bundles are the canonical models for symplectic mani- 
folds. For multisymplectic manifolds we have more variety of models (see pj) but we are 
mainly interested in the following ones. 

Let A k E be the total space of the bundle of fc-forms on a manifold E. We define a 
fc-form on A k E by setting 

9{ot)(X u ...,X k ) = a(Tp k E (a)X u Tp h E (a)X k ), 
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for a £ (A k E) x and Xx,...,X k £ T a (A k E), where p k E : A k E — ► E is the canonical 
projection. 

The form SI = —dO is a multisymplectic (k + l)-form on A k E. Of course, for k = 1 
we obtain the canonical symplectic form on A 1 ^ = T*E. 

More interesting models can be defined if E is fibered over a manifold N, say ir : 
E — ► A/ is a fibration. Indeed, we consider the bundle of r-semibasic forms on E: 

A k E = {a £ A k E | i Vl A—Av r a = 0, for r tt — vertical tangent vectors Vi,... ,v r }. 

We can reproduce the above definition to obtain a fc-form <d r on A k E which is nothing 
but the restriction of to the submanifold A k E. Therefore, the (k+ l)-form Q r — —d<d r 
is again multisymplectic and, of course, the restriction of il. 

3. Multisymplectic formalism for classical field theories 
3.1. The finite dimensional setting 

The mathematical ingredients of a classical field theory are the following ones: 

• a fibration ttxy '■ Y — ► X, where X represents a rt-dimensional space-time, that 
is, an oriented manifold of dimension n with volume form 77; 

• a lagrangian form A : Z — ► A n X, where Z denotes the 1-jet prolongation J 1 itxy- 
There are two induced fibrations ttxz ■ Z — ► X and ttxy ■ Y — > X such that 

KXZ = KXY KYZ- 

We have A = Lt], with the obvious identifications, and L : Z — ► R is the lagrangian 
function. 

In this paper we will choose fibered coordinates (x^,y l , zVj on Z such that 77 = d n x = 
dx 1 A • • • A dx n . A useful notation will be d n ~ x x l = i a n. 

Therefore, the lagrangian function is L = L{x^,y l , ). L is said to be regular if the 
hessian matrix 




We define the action integral 

M°) = [ (A)*a 

Ju 

where a is a section of ttxy defined on an open set U, and j x a denotes its first jet 
prolongation. A section a is called an extremal for the above action if 

^rM^t °<j)| t =o = 
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for every flow <p t on Y such that ttxy °<fit = Txy, and 0t(2/) = 2/ f° r au 2/ i n the boundary 
of o~(U). Since such a flow <p t is generated by a 7rjfy-vertical vector field £y vanishing on 
the boundary of o~(U), we then conclude that a is an extremal if and only if 

/ (jV)*£ Sz A = 0, (1) 
Ju 

for all satisfying the above conditions, where £z is the natural prolongation of £y to 
Z. Let us recall that if 

then 

tz = e — + . 

dy % dxt 1 dz^' 

Therefore, we deduce that a(x tM ) = (x fJ ',a l (x)) is an extremal if and only if 

Wdx^^J ? dxn = 

for all values of £*. Thus, a will be an extremal if and only if 

dL__d_dL 

dy i dx» K dzy y ' 

along j 1 a. Equations (|^) are called the Euler-Lagrange equations for L. 

Next, we will recall the multisymplectic formulation of classical field theories. 

First of all, using the volume form 77 we can construct a tensor field of type (1, n) on 

Z: 

S n = W - z^dx") A V <g> -^j. 

oz v 

The Poincare-Cartan n and (n + l)-forms are defined as 

6 L = A + s: 1 (dL), n L = -de L , 

where S* is the adjoint operator. 

Proposition 1. The (n+l)-form f2z, is multisymplectic if and only if the lagrangian 
L is regular. 

The integral action reads now as 

A{a) = / (iV)*6 i . 



The use of the Poincare-Cartan (n + l)-form allows us to characterize the extremals 
as follows. 
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Proposition 2. W a is an extremal for the action if and only if 

(j l a)*(k z Q L ) = (3) 

for all vector fields £z on Z . 

Eq. (||) gives the Geld equations, but we can consider a more general problem: to look 
for a section r of the fibration ttxz '■ Z — > X such that 

t*^ z Q l ) = (4) 

for all vector fields £z on Z. Eq. (||) is referred as the De Donder equations. 

Looking the Dc Donder problem, one realizes that we search for sections r of ttxz 
satisfying Eq. (||). We will linearize the problem as follows. First, notice that the 
subspace H T ^ = dr(x)(T x X) C T T r x \Z is horizontal, for every x in the domain of r. 
Indeed, we have 

T T ( X )Z = H T ( X ) © {Vttxz)t(x), 

where Vir denotes the vertical bundle. Therefore, these horizontal subspaces can be 
considered as infinitesimal approximations to the section in the same way that tangent 
vectors are the linear approximations to the trajectories in particle mechanics. 

Indeed, consider an arbitrary Ehresmann connection T in the fibration ttxz ■ Z — > 
X, that is, we have an horizontal distribution H which is complementary of the vertical 
one: 

TZ = H ®Vnxz. 
Denote by h the horizontal projector h : TZ — > H. We then have 

d d d d 

h(— ) = — — +r — + r — 

A direct computation shows that 



i h n L = (n-l)n L (5) 

if and only if 

(Tl -4)1 — ) = (6) 

dy* dx^ »dyid4 ^ 8^84 K " v) dy*dd { > 

(see H). 

If the lagrangian L is regular, then Eq. @ implies that r* = z^ and therefore (0) 
becomes 

8yi dx^dzl »8yi84 d^dx* ' 1 ' 
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Now, if r(x A1 ) = t 1 (x) , t^(x)) is an integral section of H we would have 

cV drt 
i _ pi f* 

" ~ dx^ ' ^~ dx v 

which proves that Eq. (||) is nothing but the Eulcr-Lagrange equations for L. 

As we have noticed, for regular lagrangians the field and De Donder equations are 
equivalent, that is, a solution r of the De Donder problem is automatically a jet prolon- 
gation of a solution of the field equations. 

Eq. (||) can be considered as the linear version of De Donder equations. If the 
lagrangian L is regular, a global solution always exists, because one can choose local 
solutions and use then a partition of the unity to glue together all the local solutions. 
But this solution will not be unique. Indeed, given a particular solution T then any other 

d 

solution is of the form h + T, where T = T l dx v ® — — is a tensor field on Z of type 
(1, n) satisfying the following conditions 

8 2 T 

T* — : r = 0, for all j. 

dzldzl 

3.2. Presymplectic formalism in the space of Cauchy data 

Next, we will show how the selection of a Cauchy surface permits to give an infinite 
dimensional setting for field theories (see || for a detailed exposition). 

A Cauchy surface is a compact oriented (n — l)-dimensional submanifold B of X. We 
shall consider a manifold of embeddings from B into X, denoted by X; it will be the 
space of parametrized Cauchy surfaces. 

A space of Cauchy data will be a manifold Z of embeddings from B into Z such that 
kxz "7 £ f° r au embedding 7 6 Z . This condition permits to define a fibration 
T xz '■ z — * x - 

A space of Dirichlet data is a manifold Y of embeddings from B to Y such that 
TT Y z 7 £ Y, for all 7 G Z. Therefore we have a fibration iTy^ : Z — > Y such that 

n xz ~ k xy n YZ- 

We define a 1-form Q on Z as follows: 

e(7)(0 - / 7*(*?e L ) 

Jb 

where 7 £ Z, and ^ G T^Z. Notice that such a tangent vector £ can be viewed as a 
mapping £ : B — ► TZ covering 7. 

Proposition 3. 0/ We have 

de( 7 )(a,6) = - / 7*(<&iei«i) 

/or a/Z tangent vectors £1, £2 G T 7 Z and /or all j G Z . 
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Let c x (t) be a curve in X. We say that c x (t) splits X if the mapping 

(b,t) ^ c x (t)(b) 

is a diffeomorphism from BxR into X. (For simplicity, we are assuming that the curves 
are defined in the whole real line). 

Let now c^(t) be a curve in Z such that its projection k X z ° c%(t) splits X. Under 
this condition, one can define a section r : X — > Z of ttxz as follows: 

t(x) = Cz(t)(b) 

since x — ■kxz ° c ^(i)(&)- 

Conversely, given a section r of 70fY and a curve cj^ (£) in X, we define a curve cg(t) 
in Z by the formula 

c^t)(b)=r(c x (t)(b)). 

Theorem 1. ^ Let t be a section of irxz, o,nd c x {t) a curve in X such that 
c x (t)(B) is contained in the domain of r. Let c^{t) be the lift of c x {t) to Z as we 
described above. If r satisfies the De Bonder equation then 

i d . {t) de = o. (9) 

Conversely, if the curve c^(t) satisfies Eq. (^) and its projection splits X, then the 
section r constructed as above satisfies the De Bonder equations. 

4. Classical field theories with non-holonomic constraints 

A non-holonomic classical field theory consists of a lagrangian function L = L(x^ , y l , z 1 ^) 
subjected to constraints of the form 

= 0, 1 < a < k. 

This means that the solutions of the field equations have to satisfy the above equations 
along the full evolution of the system. 

Definition 2. A section a of ttxy '■ Y — > A is a solution of the nonholonomic 
problem if (Q) holds for all vector fields £y on Y satisfying the condition 

<9$ Q ■ 

vrf=0, a=l,...,k. (10) 

Remark 1. The above definition could be considered as a sort of dAlembert prin- 
ciple for constrained field theories. In principle, we consider non-linear constraints. At 
the end of this section we will treat the particular case of linear or affine constraints. 
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Proceeding as in the case of unconstrained systems we deduce that a is a solution of 
the nonholonomic problem if and only if 



I 



dL _ d dL. 
dy i dx^ ^ dz l ' 



C dx n = 



for all values of satisfying (|T^) (here U is the domain of definition of a) . 

Therefore, such a solution a(x fl ) = (x^,a l (x)) would satisfy the following set of 
equations: 

dt dx^dzy dzi ( ' 

$C V>^) = o (12) 

along j 1 a. Equations (^lj) are called the non-holonomic field equations for L. Here, the 
X afx are some Lagrange multipliers to be obtained. 

We will rewrite the non-holonomic problem in a geometrical manner. 

First of all, the constraints geometrically define a submanifold M of Z. Therefore, 
a field theory with non-holonomic constraints is described by a lagrangian form A = Lt] 
and a submanifold M of Z. 

The non-holonomic De Donder equations can be then written as follows: 

i h n L - (n- 1)0 L e y(S*(TM°)), ImhcTAf (13) 

where h is the horizontal projector of a Ehresmann connection T in ttxz '■ Z — > X, TM° 
is the annihilator of TM, and y(S*(TM°)) denotes the ideal generated by the space of 
n-forms S*{TM°). 

We will check that, in fact, the above equations are correct. Indeed, the first equation 
in ( |l3| ) can be written as 

t h n L -(n-i)n L = \ a As;(dt> a ), (u) 

with some 1-forms A Q . 

In the following, we will assume the following condition: 

Assumption on the regularity of the constraints: 

The n-forms 

{S;(d$ a )\a = l,...,k} 

are linearly independent. 

Under this regularity assumption, it is esay to see that 
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Thus, if we introduce the notation 

= \ au s;(d$ a ) 

we deduce that ( |l4]) can be written as follows: 

i h fl L - (n- 1)Q L = dx 11 A 6» M . (15) 



Now, from ( |l3| ) we obtain 




(16) 



If the lagrangian L is regular, we deduce from ( |T6| ) that Tj^ — zj^, and consequently, 
( |T^ ) becomes 

If r(x M ) = (x M , t 1 (x), (x)) is an integral section of H we would have 

Br i Or 1 

» dx» ' dx v 

and thus Eq. ( O ) gives the non-holonomic field equations for L. 

Next we will give the infinite dimensional setting for non-holonomic field theory. 
To do that, we prove the following result. 

Proposition 4. Let h be the horizontal projector of an Ehresmann connection in 
the fibration ttxz ■ Z — ► X, and le.tr : X — ► Z be an integral section of this connection. 
If h is a solution of the non-holonomic De Donder equations, then we have 

T*(i ez (n L -dx»A6 tl ))=0, (19) 

for all vector fields £z on Z. Conversely, if for all the integral sections r of h Eq. ([I 
holds then h is a solution of the non-holonomic De Donder equations. 

Eq. ([l9]) can be viewed as the non-holonomic version of the De Donder problem. 

Choose now a Cauchy surface B and the corresponding spaces X, Y and Z. It should 
be noticed that now the space of Cauchy data Z consists of embeddings 7 : B — > Z 
taking values in M . 

As above we integrate the Poincare-Cartan (n + l)-form VIl to construct a 2-form 
dO. We will proceed in the same way with the constraints. Let S be the (n + l)-form 
on Z defined by 

5 = dx 1 " A0 U . 
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We then define a 2-form 3 on Z as follows: 

H(7)(£i,6) = - / 7*(<&<frH), 

for all 7 6 Z and for all tangent vectors £i, £2 6 T 7 Z- 

Theorem 2. Let t be a section of ttxz taking values in M, and c x (t) a curve in 
X such that c x (t)(B) is contained in the domain of t. Let c^{t) be the lift of c x (t) to 
Z as we described in the preceding section. If r satisfies the non-holonomic De Donder 
equation fcl<\ ) then 

i & {t) (de - a) = 0. (20) 



Conversely, if the curve cg(t) satisfies Eq. \2(\ j and its projection splits X, then the 
section t constructed as above satisfies the non-holonomic De Donder equations. 

Proof: The proof follows the same lines that in Q for the unconstrained case. Indeed, 
if r satisfies (|T^), then we have 



Cz(i)*(^2*c g (t) s ) 



= 3(cs(t))(6z(*)>6) 

= (^(t) s )(6) 

which implies (pp|). 

The converse is proved by reversing the argument. ■ 

5. Afflne and linear constraints 

Let £ be a subspace of A£ Y. As we know, an element z = j\a G Z can be viewed as 
an injective linear mapping z — da(x) : T X X — > T a r x \Y. Therefore, we could define a 
submanifold M of Z as follows: 

M = {z G Z I da(x)(T x X) is annihilated by £}. 

We will prove that M is the zero set of a family of affine constraints. 

Indeed, take a basis {(p a ,a = 1, . . . , k} of £ , and notice that each element ip a is 
locally expressed as 

tp<* = (<p a ) d n x + (ip a Y t df A d™" V. 

Since 
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where z = (x 1 - 1 , y l , z % ), we have 

<p a (da(x)-^L, . . . , da(x)^) = (<p a ) a + (^z% 

which is an afhne constraint. Linear constraints are obtained considering a subspace of 

* A?F 

the quotient vector bundle Z ' = . 

A 1 Y 

6. An example: a pneumatic tire 

We shall briefly describe a mathematical model for an automobile tire rolling with 
constant velocity, while deforming. 

There are several mathematical models: Rocard's theory, Greidanu's theory and 
Keldys's theory (see [E8|, Ch. VI, pp. 308 and ff.) The last one is the most interesting. 
We shall describe it. 

For a pneumatic tire, we introduce the following parameters describing a configura- 
tion: the coordinate x of the point K of intersection of the diameter of greatest slope of 
the wheel and the plane Oxy of the road, the angle k between the normal to the road 
and the central plane of the wheel, the angle between the Oy axis and the path of the 
central plane of the wheel on the road, the lateral displacement £ of the center of the 
area of contact of the tire with respect to K, and the angle ip measuring the torsional 
deformation of the tire (see Figure). 

The lagrangian is 

L(x,K,6,£„(p,x,k,9,£,tf) = T(x,K,e,x,k,e) -U(K,^,ip), 

where T is the kinetic energy and U is an elastic potential. The 
kinetic energy of deformation of the tire is neglected, and the 
potential energy is 

U(£, <t>, k) = i (< 2 + b(j) 2 + P Nn 2 + 2aN^n) , 

where a, 6, a, and p are constants to be experimentally deter- 
mined, and N is the normal load. It should be noticed that 
when TV increases the area of contact of the tire, which remains 
at rest during the loading of the wheel, increases. 
The kinematic constraints are 

x + i + ve + v<f) = o 

6 + <p - aV£ + 13V ip + jVk = 0, 

where a, (3 and 7 are constants, and V is the constant velocity 
of the wheel. 

The stability of vehicles on wheels with pneumatic tires is carefully analyzed in [^8| , 
introducing several simplifications to treat with the motion equations. 
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It should also be noticed that for accelerated motions, the area of the contact of 
the tire with road is not more constant (or even simply connected). This would lead 
to introduce unilateral constraints. There are mathematical descriptions for particle 
mechanic systems subjected to unilateral constraints (see [ll], || and the references 
therein) but the extension to classical field theories is not developed yet. 

7. Concluding remarks and future work 

We have developed a multisymplectic setting for first-order field theories subjected to 
non-holonomic constraints which is very similar to that for particle mechanics (see p4|). 
The dynamics in the corresponding space of Cauchy data was also analyzed. 

7.1. Algebraic computation of the Lagrange multipliers 

In Particle Mechanics the Lagrange multipliers are obtained using a simple algebraic 
method, indeed, the problem is reduced to apply the Cramer rule for systems of linear 
equations p4| . Therefore, one can use Mathematica® to do the computation. A similar 
procedure can be used for field theories, but the situation is really more complicated. In 
addition, we are trying to obtain the constrained Ehresmann connection by projecting a 
solution of the unconstrained problem. 

7.2. Applications to Continuum Mechanics 

As we said in the Introduction, the most interesting situations occur in Continuun 
Mechanics. The case of non-holonomic constraints given by connections in media with 
microstructure are being analyzed in Q| . Some independent work can be found in @@- 

7.3. Non-holonomic bracket, non-holonomic momentum mapping and re- 
duction 

These topics will be analyzed in fl and (2f| , following the lines in j7| and H . 

Acknowledgements 

The authors wish to thank the referee for his comments and remarks that help to 
improve this paper. 

REFERENCES 

[1] E. Binz, M. de Leon, D. Martm de Diego and D. Socolescu: Multisymplectic formalism for 
classical field theories subjected to non-holonomic constraints, in preparation. 

[2] E. Binz, M. de Leon and D. Socolescu: C.R. Acad. Sc. Pans, 326, Serie II b, 227-232 
(1998). 

[3] E. Binz, M. de Leon and D. Socolescu: Extracta Mathematicae 14 2, 99-125 (1999). 

[4] E. Binz, J. Sniatycki and H. Fischer: Geometry of classical fields, North-Holland, Amster- 
dam, 1988. 



[Author and title] 



13 



B. Brogliato: Non-smooth impact dynamics, Springer- Verlag, Berlin, 1996. 

F. Cantrijn, A. Ibort and M. de Leon: J. Austral. Math. Soc. (Series A) 66, 303-330 (1999). 

F. Cantrijn, A. Ibort and M. de Leon: Rend. Sem. Mat. Univ. Pol. Torino, 54, 3, 225-236 
(1996). 

F. Cantrijn, M. de Leon, J.C. Marrero and D. Martin de Diego: J. Math. Phys. 40 , 795-820 
(1999). 

G. Capriz: Continua with Micro structure, Springer Tracts in Natural Philosophy, 35, 
Springer, Berlin (1989). 

A. Ibort, E. Lacomba, M. de Leon, J.C. Marrero, D.M. de Diego and P. Pitanga: J. Phys. 
A: Math. Gen. 34, 1691-1712 (2001). 

J. Cortes, M. de Leon, D. Martin de Diego and S. Martinez: Royal Society of London 
Proceedings Series A: Math. Phys. Eng. Sci. 457 2007, 651-670 (2001). 

A. Echevarrfa-Enrfquez, M.C. Munoz-Lecanda and N. Roman-Roy: J. Math. Phys. 39, 
4578-4603 (1998). 

A. Echevarrfa-Enrfquez, M.C. Munoz-Lecanda and N. Roman-Roy: J. Phys. A 32, 8461- 
8484 (1999). 

A. Fernandez, P.L. Garcia and C. Rodrigo: Lagrangian reduction and constrained vari- 
ational calculus , in Proceedings of the IX Fall Workshop of Geometry and Physics, X. 
Gracia, M. Munoz-Lecanda and N. Roman-Roy eds., Publicaciones de la RSME vol. 3, 
Madrid, 2001. 

P.L. Garcia and A. Perez-Rendon: Comm. Math. Phys., 13 24-44 (1969). 
P.L. Garcia and A. Perez-Rendon: Arch. Rat. Mech. Anal. 43, 101-124 (1971). 

H. Goldschmidt and Sh. Sternberg: Ann. Inst. Fourier (Grenoble 23, 203-267 (1975). 

M.J. Gotay, J. Isenberg and J.E. Marsden: Momentun maps and classical relativistic fields. 
Part I: Covariant field theory, MSRI preprint (1997). 

M.J. Gotay, J. Isenberg and J.E. Marsden: Momentun maps and classical relativistic fields. 
Part II: Canonical analysis of field theories, MSRI preprint (1999). 

S.P.Hrabak: On a Multisymplectic Formulation of the Classical BRS T symmetry for First 



Order Field Theories Part I: Algebraic Structures, aath-ph/9901012 (1999) 



S.P.Hrabak: On a Multisymplectic Formulation of th e Classical BRST Symmetry for First 



Order Field Theories Part II: Geometric Structures, math-ph/9901013 (1999). 

M. de Leon: Media with micro structure: constitutive theory and dynamics, in Proceedings 
of the IX Fall Workshop of Geometry and Physics, X. Gracia, M. Munoz-Lecanda and N. 
Roman-Roy eds., Publicaciones de la RSME vol. 3, Madrid, 2001. 

M. de Leon, J.C. Marrero and J. Marin: A geometrical approach to classical field theories: 
a constraint algorithm for singular theories, in New Developments in Differential Geometry 
(Debrecen, 1994), Math. Appl. 350, Kluwer, Dordrecht, 1996, pp. 291-312. 

M. de Leon and D. Martin de Diego: J. Math. Phys. 37, 3389-3414 (1996). 

M. de Leon, D. Martm de Diego and A. Santamaria: Symmetries and reduction in classical 
field theories, in preparation. 



[Author and title] 



14 



[26] J.E. Marsden and Th.J.R. Hughes: Mathematical foundations of elasticity, Dover Publica- 
tions, Inc., New York, (1994). 

[27] J.E. Marsden, S. Pekarsky, S. Shkoller and M. West: J. Geom. Phys. 38, 253-284 (2001). 

[28] Ju.I. Neimark, N.A. Fufaev: Dynamics of Nonholonomic systems, Translations of Mathe- 
matical Monoghraphs, 33, AMS, Providence, R.I, (1972). 

[29] C. Paufler and H. Roemer: De Donder-Weyl Equations and Multisymplectic Geometry, to 
appear in Rep. Math. Phys. (2001). 

[30] S. Vignolo and D. Bruno: Iper-ideal kinetic constraints in Continuum Mechanics, preprint 
2001. 



